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We simulate quasistatic flows of an ideal two-dimensional monodisperse foam around different 
obstacles, both symmetric and asymmetric, in a channel. We record both pressure and network 
contributions to the drag and lift forces, and study them as a function of obstacle geometry. We 
show that the drag force increases linearly with the cross section of an obstacles. The lift on an 
asymmetric aerofoil-like shape is negative and increases with its arc length, mainly due to the 
pressure contribution. 
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I. INTRODUCTION 

Foams are used widely, for example in industries associated with mining, oil recovery and personal care products 
[H . Their use is often preferred because of properties such as a high surface area, low density and a yield stress 0, Q . 
In addition to this evidence of plasticity, a foam's rheology is dominated by elasticity at low strains and viscous flow 
at high strain-rates: they are elasto-visco-plastic fluids [J. 

A common probe of foam rheology is a variation of Stokes' experiment [H in which an object moves relative to 
a foam [6l-[20j. Foams have an advantage over many complex fluids in that their local structure (the bubbles) is 
observable, thus making them an excellent choice to determine the mechanisms by which non-Newtonian fluids show 
different responses to Newtonian fluids. In addition, a two-dimensional foam is a realizable entity, for example the 
Bragg bubble raft [2]| , with which it is possible to perform a rheological experiment in which the shape and velocity 
of each bubble can be tracked in time. Foams are also amenable to numerical simulation because of the precise local 
geometry that is found wherever soap films meet. Plateau's laws, which describe how the films meet, are a consequence 
of each soap film minimizing its energy, equivalent to surface area, and it is this that provides the algorithm for the 
work described here. 

For flow to occur in a foam, the bubbles must slide past each other. This occurs through Tl neighbour switching 
topological changes [l^l, in which small faces and/or short films disappear and new ones appear. Sometimes referred 
to as plastic events, these are a visible indication of plasticity in a foam, and act to reduce the stress and energy. 
Numerous contributions to viscous dissipation occur 2^ , although we assume that if the fiow is slow enough they 
can all be neglected. 

DoUet et al. [l5| measured the drag, lift and torque on an ellipse in a two-dimensional foam flow in a channel. The 
lift was maximized when the ellipse was oriented at an angle of 7r/4 to the direction of flow. DoUet et al. Q found 
that an aerofoil embedded in a foam flow exhibited a negative lift, which they attributed to the elasticity of the foam. 
This augments the list of well-known non- Newtonian effects that contradict the sense of what is known for Newtonian 
fluids. 

We present here elasto-plastic simulations, in the so-called quasi-static limit, for 2D foam flow around an obstacle, 
and investigate the effect of the symmetry of the obstacle in determining the magnitude and direction of the drag 
and lift. Such simulations allow us to exclude consideration of viscous effects, and even to separate out pressure and 
fllm network contributions to the forces on an obstacle, both of which are difficult to do in experiment. As a means 
of determining drag and lift on an obstacle, they have been validated against experiments on an ellipse [l5| by Davies 
and Cox H. 

We consider a range of obstacle shapes, illustrated in figure[TJ Since the Evolver uses a gradient descent method, we 
are unable to simulate an obstacle with sharp corners. We therefore round the corners of each obstacle with segments 
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FIG. 1: Pictures of the obstacles, oriented with flow from left to right: (a) circle, (b) horizontal stadium, (c) square, (d) 
symmetric aerofoil and (e) asymmetric aerofoil. 



of a circle to smooth the boundary. The shapes are: 

(a) a circle, which provides the standard case with full symmetry. Its cross-section is H = 2R. 

(b) the union of a square and two semi-circles, which we call a "stadium" , arranged either vertically or horizontally. 
The side-length 2R of the square is equal to the diameter of each semi-circle, so that the area is determined by 
just one parameter, R. The cross-section is 2R (horizontal stadium) or 4i? (vertical stadium). 

(c) a square, with rounded corners. The radius of curvature of the corners is set to one-eighth of the side-length of 
the square, R — L/8 , so that the area is again determined by just one parameter, and H — L. Also a diamond, 
which is the square rotated by 7r/4, with H sa ^/2L. 

(d) a symmetric aerofoil, with long axis parallel to the direction of foam flow, defined by two arcs of circles bounded 
by two tangential straight lines. Three parameters are needed: length L (distance between the centres of the 
circles), and radii Ri (leading edge) and R2 (trailing edge). This shape has up-down symmetry but not fore-aft 
symmetry, and cross-section H = 2max{Ri, R2). If Ri = R2, then this is a "long" horizontal stadium. 

(e) an aerofoil-like shape with up-down asymmetry, in which two circles of equal radius R2 are joined by arcs of 
radius i?i and Ri + 2R2. The distance between the circles is parametrized by the angle 9i. Its cross-section is 
H = {Ri + R2){1 — cos 6*1) -I- 2i?2- This approximation to a standard aerofoil dispenses with the singular point at 
the trailing edge. 

We begin by describing our numerical method (211) • The forces on each obstacle are given in ^III At we find that 
the drag is mainly determined by its maximum cross-section H perpendicular to the direction of flow and that a 
significant lift is found only for the aerofoil without up-down symmetry. The field of bubble pressure around the 
obstacle, which is the main contribution to this lift, is described in ijlll B[ and we make some concluding remarks in 
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FIG. 2: (a) Sketch of the simulation, in tliis case for a horizontal stadium. A 2D foam is created between two fixed walls 
and caused to flow in the positive x direction by increasing the area of the region to the left of the dark line of films joining 
the two walls. The obstacle is created in the centre of the channel; each film that touches the obstacle applies an equal force 
outward in the direction normal to the obstacle and each bubble applies a pressure force inward at the middle of the shared 
boundary. The films bunch up at the trailing edge of the obstacle and the bubble pressures rise at the leading edge due to the 
flow, leading to drag and lift forces on the obstacle, (b) Example (vertical stadium, area ratio = 6) of the pressure [Fp) and 
network (-Ft) contributions to the drag (a;) and the lift (j/) as a function of iteration number. The drag forces increase linearly 
before developing a saw-tooth variation which is linked to a build-up of stress followed by avalanches of Tls in the foam. The 
horizontal lines show the average drag forces. In this case the pressure and network contributions to the lift are both negligible. 



II. METHOD 

We use the Surface Evolver [2^ in the manner described by Davies and Cox [2^ . We create three foams of around 
725 bubbles (in this range the number of bubbles does not affect the results; data not shown) between parallel walls 
with a Voronoi construction [l^[2^. The channel has unit length and width W = 0.8. The foams are monodisperse, 
with bubble area denoted Ai, and about 22 bubbles in the cross-section of the channel. A bubble in the centre of 
the channel is chosen to represent the obstacle, and its periphery constrained to the required shape; its area is then 
increased until it reaches the desired area ratio = Aobs/Ai, and it is then fixed - see figure [2Ja). The tension of 
each film, 7, which is twice the air-liquid surface tension and is in effect a line tension, is taken equal to one, without 
loss of generality. 

The boundary conditions are that of free slip on the boundary of the obstacle and the channel walls, so that the 
films meet the boundaries at 90°, and periodicity in the direction of flow. We checked in a few instances that changing 
the boundary condition on the channel walls to non-slip has little effect on the forces on a small obstacle in the centre 
of the channel. At each iteration the foam is pushed with a small area increment dA = 5 x 10~^ to create a pressure 
gradient p^ . The perimeter is then evolved towards a local minimum and Tls are performed whenever a film length 
shrinks below = 1 x 10"'^ (representing a foam with low liquid fraction, of the order of lO"**). A simulation runs for 
1500 iterations to ensure that the measurements are made beyond any transient in which the foam retains a memory 
of its initial state. Each simulation takes about one week on a 1.5GHz CPU. The method has been validated against 
experiment in the case of an elliptical obstacle [l^ [IJ] . 



A. Drag and lift 

Each film that touches the obstacle applies an outward force with magnitude equal to the force of surface tension 
and direction perpendicular to the obstacle boundary. Their resultant is the network force 

FT^iY^n, (1) 

where Hi is the unit outward normal at the vertex i terminating each film that meets the obstacle. See figure [21[a). 
Each bubble that touches the obstacle applies a pressure force inward at the middle of the shared boundary. Their 
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resultant is the pressure force 

Fp^-Y^Pjl^rij (2) 

j 

where pj is the pressure of bubble j, Ij the length of shared boundary and fij the unit outward normal to the obstacle 
at the midpoint of the line joining the two ends of the shared boundary. 

The drag on an obstacle is the component of the sum of the network and pressure forces in the direction of motion, 
Fd — + Fp. The lift is the component perpendicular to this, Fl = Fl^ + Fp, with the convention that positive 
values of lift act in the positive y direction. All four components are recorded at the end of each iteration, and averaged 
above 600 iterations, well beyond any transient. An example is shown in figure [^^b). The standard deviation of the 
fluctuations in force about this average are used to give the error bars in the figures below. 

III. RESULTS 
A. Drag and lift force on an obstacle 

The drag and lift oscillate in a saw-tooth fashion (figure [2jb)), caused by intervals in which the imposed strain is 
stored elastically followed by cascades of Tl topological changes. Nonetheless, they have a well-defined average. We 
find that for all obstacles with up-down symmetry the average lift is close to zero. 

We vary the area ratio of each obstacle, usually in the range one to ten but o ccasiona lly higher. We normalize the 
cross-section and length of each obstacle by the average bubble diameter db = \f^A\^ which, since the walls are far 
enough away not to have an effect on the drag and lift, is the significant length-scale here. We choose to plot the 
resulting drag as a function of cross-section Hj db (figure [3]) since it gives an approximately linear relationship . It 
is apparent that the drag increases with obstacle cross-section most quickly for "blunt" objects with a vertical leading 
edge (square, vertical stadium). Obstacles with a rounded leading edge (circle [l^l, horizontal stadium) experience 
lower drag for given cross-section. In each case, the main contribution to the drag is usually due to network forces; 
the pressure contribution to the total drag is lower but follows the same trends. 

To tease out the effect of obstacle shape on the two components of drag studied here, we fix the cross-section (figure 
|3]Ja),(b)) and vary the shape. The pressure contribution to the drag is highest when the leading edge is blunt (vertical 
stadium, square), since this causes the greatest deformation to the bubbles. Similarly, the network contribution to 
the drag is highest when the trailing edge is rounded (the most "circular" case in figure |U[b)), although this effect is 
weaker, since a rounded trailing edge allows more films to collect in that area. The shape of the diamond is such that 
the network drag is very low, since films can gather on the sloping sides as well as the rounded region at the very tip 
of the trailing edge, while the pressure drag is intermediate. 

The length L of an obstacle has only a weak effect on the drag (figure |U[c)). In particular, this is the case for a 
symmetric aerofoil with i?2 = Ri, since most of the films that touch the obstacle are perpendicular to the direction 
of foam flow. By varying the ratio R2/R1 for a symmetric aerofoil with fixed cross-section H and fixed area ratio 
Qr = 10, we can investigate the effect of fore-aft asymmetry. Figure |l{d) shows that the total drag varies little, 
emphasizing that cross-section and rounded leading and trailing edges make the major contribution to the drag. The 
pressure contribution to the drag decreases with R2/R1, that is, as the leading edge gets smaller and bubbles are less 
deformed there. 

The lift is, on average, zero for all obstacles with a horizontal axis of symmetry (as in figure O^b)); it is only 
significant for the asymmetric aerofoil, being negative and of the same order of magnitude as the drag. In particular 
the lift increases with aerofoil length (figure [5]), and the major component of lift arises from the bubble pressures. It 
appears therefore that the curvature of the aerofoil induces changes in bubble pressures, and that it is this, rather 
than an imbalance in the number of films pulling on the top and bottom surfaces of the object, that gives rise to the 
lift. We return to the bubble pressures below. 

To test the effect of obstacle position in the channel, we placed the same asymmetric aerofoil in three different 
positions across the channel: y = 0.251^, 0.514^ (reference case) and 0.751^. No significant difference in the drag or 
lift was observed (data not shown), indicating that the obstacle was still sufficiently far from the walls that they don't 
interfere with the flow (recall that this is a elasto-plastic rather than a viscous flow, distinct from a Newtonian fluid 
where the wall always has an effect in 2D) and that the lift is not just due to the foam squeezing through the gap 
between wall and obstacle. 
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FIG. 3: (Color online) Drag vs obstacle cross-section H/dt. Images are for obstacles with area ratio = 10 with flow from 
left to right, (a) Vertical stadium (a^ = 2,3,4,6,8,10). (b) Horizontal stadium (a,. = 2,4,5,6,8,10,20,30). (c) Square 
(flr = 2, 4, 6, 8, 10, L/R = 8). (d) Diamond (a^ = 2, 4, 6, 8, 10, L/R = 8). 

B. Pressure field around an obstacle 

To further probe the phenomenon of negative lift in foams, in figure [S] we compare the distribution of bubble 
pressures around the flat-bottomed aerofoil with an up-down symmetric obstacle typified by the square. The Surface 
Evolver calculates the bubble pressures (as Lagrange multipliers of the area constraints) in such a way that they are 
all relative to the pressure of one bubble. Thus the average pressure is subtracted from all values at each iteration, 
before binning the data as above. 

The bubble pressures decrease in the x direction, on average, because of the flow. The presence of an obstacle 
induces a region of high pressure at the leading edge and a region of low pressure at the trailing edge. In addition, 
the asymmetric aerofoil shows a region of high pressure above and low pressure below, confirming that the pressure 
contribution to the lift is downwards. 



IV. CONCLUSIONS 

The simulations described here show that the forces on an obstacle embedded in a flow of foam depend strongly 
on the shape of the obstacle. We separate two components, due to the pressure in the bubbles and the network of 
soap films, and find that the pressure contribution decreases with the rounding of the leading edge and the network 
contribution decreases with the rounding of the trailing edge. Further evidence is given in figure [Tl^a). 

In classical fluid mechanics, the presence of viscosity can give rise to trailing vortices and circulation around an 
obstacle in a fluid flow. Here, not only do we neglect viscosity, but the discrete nature of the foam probably suppresses 
any possibility of circulation. Yet a lift force is still observed for obstacles without lateral symmetry, and it arises 
because of the way in which the obstacle deforms the bubbles that make up the foam. It is therefore an effect of 
elasticity or, more generally, viscoelasticity [13, HI], due to the normal stresses generated in the fluid, and acts in 
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FIG. 4: Drag force on different obstacles, (a) Drag vs shape at constant cross-section H/db ~ 2.1. The pressure contribution 
to the drag decreases with the rounding of the leading edge and the network contribution decreases with the rounding of the 
trailing edge, (b) Drag vs roundness R/{L + R), interpolating between a square (7i = L/8) and a circle {L = 0) with Ur ~ 10. 
The same effect is seen as in (a), (c) Drag vs obstacle length, measured as (L + 2R)/di,, for symmetric aerofoils with Tii = R2 
at constant cross-section H/di, « 2.1. The first point on the left corresponds to a circle (L — 0), and the second to a horizontal 
stadium (L — 2R). The network contribution to the drag decreases slightly with length, (d) Drag vs radius ratio R2/R1 for a 
symmetric aerofoil (a^ = 10, L varies). The pressure drag decreases when the leading edge has a smaller radius of curvature. 



the opposite direction to the usual sense of "lift" . A concave underside, as in the famihar Joukowski profile and the 
asymmetric aerofoil described above, is not necessary to obtain a negative lift (figure [7]Jb)). 

It remains to determine whether a given obstacle is actually stable with respect to rotation; that is, whether the 
torque on any given obstacle is sufficient to rotate it and thereby reduce the drag and/or lift. This is a necessary 
pre-cursor to using this work to determining which shapes of obstacles offer the least resistance to foam flow. It is also 
of interest to incorporate some element of viscous dissipation, perhaps using the viscous froth model [2^ , within the 
simulations, which has a particularly significant effect on rotation [24| but also the film motion around an obstacle. 
We shall return to both these issues in future work. 
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(c) 

FIG. 6; Pressure fields averaged over the duration of the simulation, (a) Square obstacle with ar = 10, showing increased 
pressure upstream of the obstacle and low pressure downstream, (b) Asymmetric aerofoil, with iJi/^/d^ = 3, Rijyfd^ = 0.75 
and Q\ — n/6, showing low pressure beneath as well as downstream. The increase of pressure upstream is less-pronounced, and 
there is a pressure peak beneath the trailing edge of the aerofoil, (c) Zoom of the typical arrangement of films around the same 
aerofoil, with bubbles shaded by instantaneous pressure on a scale by which pressure increases with grey intensity. In both 
representations a region of low pressure is evident beneath the aerofoil - it is this which induces a negative lift - as well as the 
pressure peak beneath the trailing edge. 
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FIG. 7: (Color online) (a) An interpolation between a diamond and a circle takes the shape shown, with = 10. When its 
sense is flipped relative to the direction of flow, the relative contributions to the pressure and network drag change, while the 
total drag remains the same: a sharp leading edge and rounded trailing edge reduces the pressure drag and increases the network 
drag, (b) Instantaneous arrangement of films around a flat-bottomed aerofoil (parameters: Ur = 9, cross-section H/dt = 1.59, 
radius of curvature of leading edge is Ri/di, = 0.41, of trailing edge is R2/db = 0.17 and of upper side is Rs/di, = 3.46). The 
instantaneous values of drag and lift are Fy = —2.90, = —2.36, F^ = 2.03, F^ = 0.42. The lift is again negative, both 
network and pressure contributions are similar, and the total lift is of the same order of magnitude as the total drag. 



